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Abstract. On the basis of detailed first-principles calculations the anisotropic 
thermoelectric transport properties of biaxially strained silicon were studied 
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powerfactor was only found for the hole-doped case under small biaxial tensile 
strain. The results are discussed in terms of band-structure effects. An analytical 
model is presented to understand the fact that the thermopower decreases if 
degenerate bands are energetically lifted due to a strain- induced redistribution 
of states. 
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1. Introduction 

Thermoelectrics convert heat into electric current, and vice versa. Known for more 
than 60 years, thermoelectrics attract currently attention [1, 2, 3, 4]. With nearly 
90 per cent of the worlds power being generated by low efficient heat engines, 
thermoelectric modules could potentially convert parts of this wasted heat into 
electricity. Their conversion efficiency can be stated by the figure of merit 

rTS'2 

ZT = ———T, (1) 

where a is the electrical conductivity, S the thermopowcr, k^i and Kph are the electronic 
and phononic contribution to the thermal conductivity, respectively. The numerator 
of Eq. 1 is called powerfactor PF = aS'^ . 

While ZT > 1 was challenging to be reached in the last decades, 
nowadays nanostructured thermoelectrics enable even larger values of ZT [5, 6, 7]. 
Unfortunately, those materials arc often based on environmentally unfriendly lead, 
tellurium or selenium compounds and arc therefore hard to integrate in semiconductor 
electronics. 

Silicon, the cradle of modern semiconductor electronics, is nonpoUuting, readily 
available, cheap and perfectly integrated in present electronics infrastructure. While 
silicon has been stated as inefficient thermoelectric in the past due to his enormous 
thermal conductivity [8] , recent experiments and theory revealed that nanostructuring 
could lead to thermoelectric efficiencies comparable to state of the art commercial 
thermoelectric materials [9, fO, 11, 12, 13]. 

As in nanostructured materials mechanical strain plays an important role, this 
paper will focus on the influence of biaxial strain on the electronic thermoelectric 
transport properties of bulk silicon , which might occur in rolled up and layered Si 
heterostructurcs [14, 15, 16, 17]. It is well known, that similar strain physics, e.g. band 
splitting and band deformation are expected for strained three dimensional systems, 
as well as for one and two dimensional silicon devices in complex strain states. A 
comprehensive overview on this can be found in ref. [14]. In detail these common 
strain physics affect the electronic carrier transport depending on the dimensionality, 
temperature and doping density. The question to be answered in this paper will be 
whether tensile or compressive strain will lead to an enhancement or reduction of the 
powerfactor in silicon under a certain doping and temperature environment. Beside 
our interest in the high temperature thermoelectric application of strained silicon we 
want to emphasize the possible importance of our results in the low temperature 
regime for the mctal-oxidc-semiconductor device community, where the knowledge of 
the thermoelectric properties of silicon under strain could help to understand parasitic 
effects in these devices. In the low doping regime at low temperature, an enhancement 
of one part of the powerfactor, namely the electrical conductivity, under externally 
applied strain was found and heavily investigated in the last decades. [18, 19, 14] For 
this purpose the paper will be organized as follows. In section 2 we introduce our 
first principle electronic structure calculations based on density functional theory and 
the transport calculations based on the solution of the linearized Boltzmann equation. 
With this knowledge we start the discussion of the thermoelectric transport properties 
of unstrained bulk silicon (Sec. 3.1) and present afterwards the influence of biaxial 
strain on the electron or hole doped case of silicon in section 3.2. In the last paragraphs 
3.3 and 4 the optimal powerfactor under strain due to variation of doping is determined 
and analyzed, while estimations of the possible figure of merit are given. 
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2. Methodology 

Our approach is based on two ingredients: first principles density functional theory 
calculations (DFT), as implemented in the QuantumEspresso package [20] and an 
in-house developed Boltzmann transport code [21] to calculate the thermoelectric 
transport properties. 

2.1. Electronic structure 

In a first step the band structure of the strained and unstrained Si was calculated 
using the general gradient approximation (GGA) and the PBE exchange correlation 
functional [22]. Fully relativistic and norm-conserving pseudo potentials [23] were 
used to treat the spin-orbit splitting of the Si valence bands in a correct way. The 
calculations were performed with the experimental lattice constant a = 5.434 A for a 
face-centered tetragonal eight atom unit cell. The in-plane biaxial strain is simulated 
by changing the ^^/a ratio but keeping the cell volume constant. Throughout the paper 
the biaxial strain will be given in units of the relative change of the in-plane lattice 
constant as ^°/ao = "/ao — 1- That means tensile strain considers changes '^"■/ao > 0, 
while compressive strain means ^"-/ao < 0. 

As expected, our DFT calculations underestimate the size of the band gap at 
zero temperature and do not reproduce the temperature dependence of the gap. For 
this purpose we included a temperature-dependent scissor operator [24], so that the 
strain- and temperature-dependent energy gap Eg in eV becomes 

Eg{T,—) = Eg{T = 0,—) + UGGA (2) 
ao ao 

4.73 • 10"'* • T2 

r-)-636 ' 

where Eg{T = 0, ^) is the zero temperature gap obtained by our self-consistent DFT 
calculations, J/qga = 0.57 eV is a static correction to fulfill the experimental low 
temperature gap. The third part of equation 2 is the correction of the temperature 
dependence of the band gap [25] in a wide temperature range with T given in K. 

2.2. Boltzmann transport 

With the converged results from the first step we are now able to obtain the 
thermoelectric transport properties by solving the linearized Boltzmann equation 
in relaxation time approximation (RTA) [21]. Boltzmann transport calculations for 
thermoelectrics have been carried out for quite a long time and show reliable results 
for metals [26, 27, 28, 29] as well as for wide- and narrow gap semiconductors [30, 31, 
32, 33]. Within in here the relaxation time is assumed to be constant with respect 
to wave vector k and energy on the scale of ksT. This assumption is widely used for 
metals and doped semiconductors. The constant relaxation time is a big advantage 
for the calculation of the thermopower S without any adjustable parameter, while 
lacking any doping or temperature dependence of t. For unstrained silicon, doping 
dependent relaxation times in the order of 15-150 fs for electron doping and 6-65 fs for 
hole doping could be estimated from experiments in ref. [34]. To concentrate on the 
bandstructure effects we assume the relaxation time does not depend on strain, while 
it was shown, that under strain the dominant scattering process varies: For unstrained 
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Si, the room-temperature scattering is dominated by optical phonons, i.e., intervalley 
scattering, whereas for strained Si, this scattering process is less efficient [35, 36]. 

The temperature- and doping- dependent thermopower in- and cross-plane is 
defined as 

1 c^l\{^i.T) 



where 



eT40)(^,r)^ 



(4) 

is the transport distribution function as termed by Mahan and Sofo [37], for given 
chemical potential fi at temperature T and carrier concentration N determined by an 
integration over the density of states n{E) 

YB"^'"' fj.+AE 

J dEn{E)[f^^^T)-l]+ J dEniE)f^^^T), (5) 

li-AE CB""'" 

where CB'"™ is the conduction band minimum and VB""^^ is the valence band 
maximum. As can be seen straight forward from eqs. 3 and 4 the electrical conductivity 
(7 is then given by 

a,_„ =2e24°)|(A.,T). (6) 

The Fermi-surface integration, which is requested in equation 4, is performed 
within an extended tetrahedron method [38, 39, 40] interpolating the calculated 
Eigenvalues E^ on a mesh of at least 13950 k points in the irreducible part of the 
Brillouin zone. ^^±\\ (^i T ^ 0) was determined on a fine energy mesh with a step 
width of 1 meV. Convergence tests confirmed that the calculation of the electrical 
conductivity and the thermopower at a given temperature via equations 6 and 3 
requires {E, T = 0) for a quite large range. Convergence of the integrals 4 and 5 
was achieved with an adaptive integration method for 2AE > 2QkbT in the limit of 
low carrier concentrations TV < 1 • lO^^'cm"'^. 



3. Thermoelectric transport 

3.1. Unstrained case 

Since the thermoelectric transport properties of the strained silicon will always be 
discussed in comparison with the unstrained case, we will provide an insight into the 
transport properties of the unstrained silicon, first. In figure 1, the thermopower and 
the powerfactor are shown for two different temperatures (500K and 900K) in a wide 
doping range. The qualitative behavior of the electrical conductivity is indicated by 
the dashed-dotted line to emphasize the trend of the resulting powerfactor. The picture 
is well known for the interrelation of electronic transport and the thermoelectric 
properties of semiconductors [41, 1]. At constant temperature the thermopower (cf. 
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Figure 1. (color online) In-plane thermopower (red dashed lines, refer to left 
scale, red units) and powerfactor (blue solid lines, refer to right scale, blue 
units) for unstrained, electron doped silicon in dependence on the doping level. 
Furthermore the doping dependence of the electrical conductivity is given as a 
dashed dotted line in arb. units. The maxima of the powerfactor are marked by 
black open circles. 



red dashed lines in figure 1) decreases at very low and higher doping level and reaches 
a maximum inbetween [41]. 

As can be seen from figure 1 the thermopower reaches a maximum around an 
electron carrier concentration of 1 • 10^^ cm^^ at 500K, while the maximum at 900K 
is shifted to a larger doping level of 1 • lO^'^ cm~^. Beside that, the maximum of the 
more relevant powerfactor (cf. blue solid lines in figure 1, optimal values indicated 
by black open circles) is shifted to huge electron carrier concentrations of about 
1 • 10^^ cm~^. This is determined by the linear increase of the electrical conductivity 
with increasing charge carrier concentration. Obviously, there is not much space to 
optimize the powerfactor with respect to temperature and charge carrier concentration 
for unstrained silicon. We will focus on this optimization in more detail in subsection 
3.3. The complex dependence of the thermopower on temperature and doping is 
shown in figure 2 for electron and hole doped unstrained silicon at different doping 
concentrations. For low temperatures and doping levels the thermopower reaches 
values of 1000 ^V/K and above, which is caused by the location of the chemical 
potential near the band edges. The denominator of equation 3, proportional to the 
electrical conductivity, is small, while the nominator is large, the thermopower gets 
maximal. At fixed charge carrier concentration the position of the chemical potential is 
shifted towards the middle of the gap with increasing temperature. The nominator in 
equation 3 decreases, because it is mainly determined by the opposite contributions of 
the tails of the derivative of the Fermi-Dirac distribution function with respect to the 
valence and conduction band (equation 4, n=l). At a distinct temperature of about 
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Figure 2. (color online) Thermopowcr for unstrained silicon in dependence 
on doping and temperature. Electron doping refers to the blue lines in the 
lower part of the figure, while red lines refer to hole doping and positive values 
of the thermopower. Experimental data (squares, circles and triangles) from 
ref. [42, 43, 44] are given for comparison. 

500 K the thermopowcr rapidly vanishes. At this temperature the electronic transport 
enters into the bipolar intrinsic regime. To emphasize the relevance of our calculations 
experimental results for merely pure silicon in the intrinsic transport regime arc added 
in figure 2. We want to point out, that the thermopower does not vanish at all, but 
converges to small negative values for electron, as well as for hole doping. At large 
charge carrier concentrations of about 0.16 e/uc (dotted lines in figure 2), where the 
powerfactor becomes large, the thermopower grows linearly with temperature upto 
values of 150/iV/K at 1000 K. In the heavy doped regime the chemical potential is 
located deeply in the bands and equation 3 qualitatively coincides with the well known 
MOTT relation S oc ^^'^/J^^^ for the thermopower in RTA [45]. 

3.2. Optimization of the powerfactor by strain 

3.2.1. Electron doping Having provided a general view on the thermoelectric 
transport properties of unstrained silicon above, we will now focus on the effect of 
biaxial strain on those properties. The results are presented in comparison to the 
unstrained case, starting with the electron doped case, followed by the hole doped 
case in subsection 3.2.2. 

In strain- free bulk silicon, as introduced in section 3.1, the conduction-band 
minimum (CBM) consists of six equivalent valleys on the F-X high symmetry line as 
shown in the middle panel of figure 3(b). The Fermi surface pockets corresponding to 
these valleys are shown on top, with the absolute value of the carriers group velocity. 
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Figure 3. (color online) Fermi surfaces of electron doped silicon under 
compressive strain (left), no strain (middle) and tensile strain (right). On the 
Fermi surfaces the absolute value of the group velocities are plotted in units of 
0.08 • 10® m/s. As reference the band structure on two high symmetry lines is 
given below. The doping corresponds to additionally 0.01 electrons per unit cell 
which causes carrier densities of 6.25 ■ 10^^ cm"''. 

Vk entering equation 4, plotted in color on the surface. The color code gives Vk in 
units of 0.08 • 10® m/s. By applying biaxial in-plane strain, the six CBM valleys 
are energetically split into two groups. Four degenerate in-plane A4 valleys and two 
degenerate cross-plane A2 valleys (see figure 3 (a),(c)). While the number of electrons 
is fixed, the different sizes of the ellipsoidal pockets is caused by a change of occupation 
numbers under strain. The color code indicates the overall smaller Fermi velocities 
on the small pockets and in particular on the principal axis of the pockets, whereas 
larger velocities are found for states propagating perpendicular to the principal pocket 
axis. In figure 4 the thermoelectric transport properties of biaxial strained silicon 
for two fixed electron doping regimes are shown. First, we consider the left column 
(figures 4(a)-(c)), which shows the electrical conductivity a, the thermopower S and 
the resulting powerfactor PF ~ aS^ for a low electron doping of 1.6 ■ 10~^e/uc and a 
low temperature of 100 K. For this doping and temperature value a large enhancement 
of up to 35% of the electrical conductivity is obtained for the in-planc component a^^ at 
low tensile strain and for the cross-plane component cr^at small compressive strain, 
cr^ drops noticeable under small tensile strain up to 30% of the unstrained case, while 
experiences a slight drop down to 83% under small compressive strain, respectively. 
With the conduction bands of silicon behaving parabolically near the band 
edges, the calculated transport properties can be understood for small charge carrier 
concentrations and low temperatures in terms of effective masses and occupation 
number redistributions (see [35]. With increasing tensile strain the A4 bands lift 
up and the occupied states from the higher bands are transferred to the lowered A2 
bands (comp. figure 3(c)). At a certain tensile strain the A4 pockets are completely 
depleted and solely the maximally occupied A2 states contribute to the transport. 
In simplified consideration one can estimate the relative change in the electronic 
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Figure 4. (color online) Anisotropic thermoelectric transport properties for fixed 
temperature and electron doping concentrations in dependence on compressive 
and tensile strain. Left panels (a)-(c)) correspond to an electron doping of 
1.6 • 10~^ e/uc at a temperature of 100 K, while right panels refer to an electron 
doping of 0.48e/uc at a temperature of 900 K. On the left axis of each figure the 
relative value compared to the unstrained case is shown, while on the right axis 
the absolute values are given. Electrical conductivity ((a),(d)) and powerfactor 
((c), (f)) are presented in units of the relaxation time r. 



conductivity from the relative change in the effective electron mass. As reported 
earlier [35] it is toj_,a2 = ™j_,A4 = 0.205 and ?ti||,A2 = ™||,A4 = 0.926, whereas the 
masses arc in units of the free electron rest mass. For the in-plane component cr^^ of the 
electrical conductivity at sufficient tensile strain only the lowered A2 bands contribute 
with their perpendicular mass m±^A2- With the notation 
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0.277' 



(7) 



(8) 



the normalized asymptotic value becomes: 

1 / 2 \ 

- • TOO = 1-35 

2 V'^-L,A2/ 

For the in-plane component a^^ of the electrical conductivity at compressive strain only 
the four pockets of the lowered A4 bands contribute equally with their parallel and 
perpendicular mass: 

if— + — 

4 V'"||,A4 TO^,A4 



■ TOO = 0.83. 



(9) 
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Figure 5. Analytical dependence of the thermopowcr on biaxial strain for 
different electron charge carrier concentrations, small charge carrier concentration 
(dashes lines) and increased charge carrier concentration by a factor five (solid 
lines). 



For the cross-plane conductivity cr^at tensile strain it is 

- mo = 0.30, (10) 

whereas at compressive strain cTj^ becomes 

U^-) - mo = 1.35. (11) 

Since the powerfactor is composed out of a and S we analyze the influence of strain 
on the thermopower as well. In the low temperature and doping regime (figure 4 left 
panel) no enhancement of thermopower at either compressive or tensile strain could 
be found. It can bee seen, that for tensile strain the thermopower decreases by about 
10%, while for compressive strain a drop of about 5% is found. The difference between 
the in-plane and cross-plane thermopower is marginal. In terms of Motts formula 
[45] it means that the energy dependence of jCf^^^ {E, T) and C^l\^ {E, T) is almost the 
same. Nevertheless it is interesting to understand why the thermopower of silicon is 
reduced under strain and why the decay changes for tensile and compressive strain. 

For qualitative understanding of our ab initio results we apply a free electron 
model to discuss the thermopowers behavior on biaxial strain. The strain dependent 
electrical conductivity at zero temperature was modeled as proposed by [35] and then 
the thermopower was calculated by Motts relation in RTA [45]. Figure 5 shows 
the resulting thermopower under tensile and compressive strain for a small (dashed 
line) and a five times larger charge carrier concentration (solid line). For small charge 
carrier concentration the thermopower rapidly drops to constant values of 48% for 
tensile and 76% for compressive strain, respectively. We note that the maximum of 
S which is increased by 4% compared to the unstrained case is not located at zero 
strain and is shifted to very small values of compressive strain. This behavior is more 
pronounced for the large charge carrier concentration, where the maximum of the 
thermopower becomes wider and is shifted to values of -0.25% compressive strain. 
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Again the enhancement reached for S is about 4% and the asymptotic values keep 
unchanged with doping. We note that for the calculation of the thermopower the 
influence of the effective masses completely cancels. For this reason neither for the 
ab initio nor the analytic calculation a noticeable difference of the in-plane and cross- 
plane thermopower was found. The different saturation values of the thermopower 
under tensile and compressive strain again however can be explained in terms of 
a redistribution of electrons. For free electrons at T=Q the electrical conductivity 
a oc E^^^ and the resulting thermopower S cx E^^. The amount of charge carriers 
in one spin band is given by iV oc E^/-^. With Motts relation S oc N~'^/^. As 
stated before in this paper under sufficient tensile strain only the twofold degenerated 
A2 pockets are occupied compared to the sixfold degenerated CBM pockets in the 
unstrained case, so the occupation of every pocket increases by a factor of 3. One 
directly yields St„„../so = (|)~^/^ = 0.48. For compressive strain only the fourfold 
degenerated A4 pockets contribute, so it is -Scompyso = (f)^^^'^ = 0.76. The fact, 
that for larger charge carrier concentration larger strain has to be applied to reach 
these limits, is linked to the fact that higher strain is required to reach a complete 
redistribution of states into either the A2 or A4 pockets. As consequence of the 
discussed results we see in figure 4(c) an enhancement of the powerfactor in cross- 
plane direction up to 27% at small compressive strain, while the in-plane powerfactor 
is only marginally enhanced by about 5% under low tensile strain. In-plane transport 
under tensile strain at low doping and low temperature plays an important role in 
silicon based devices. Within figure 4(c) it is obvious that the strain induced influence 
of the powerfactor on this transport will play a minor role. We want to point out that 
the results on the thermopower discussed above arc generally valid for all systems with 
degenerate occupied bands. Lifting of the degeneracy causes redistribution of electrons 
and reduction of the thermopower. While the low temperature and low doping case was 
convenient to provide some general findings on an analytical level, we will now focus on 
the high temperature and high doping regime (see figure 4 left panels) where the power 
factor might be suitable for thermoelectric application (see also 1). At a temperature 
of 900K the electronic bandstructure on a width of at least ±A£' = ±770 meV around 
the position of the chemical potential has to be included, which makes a description of 
the electronic transport properties within a spherical band picture impossible. Rather 
than providing analytical quantities a more qualitative description of our ah initio 
calculations will be given instead. The electrical conductivity in figure4(d) states the 
same qualitative tendencies for (T||and a^as derived for the low temperature case. As 
a consequence of the high temperature and the related broadening of the Fermi-Dirac 
distribution in equations 4 and 5 as well as the high charge carrier concentration the 
redistribution of states as described before is not completed for the strain values 
considered here. The analytical limits for the enhancement of errand cr^for the 
given high doping and temperature should be achieved for biaxial strains of at least 
^a/ao = ±3%, respectively. It is worth mentioning, that the absolute value (cf. right 
scales in figures 4(a), (d)) of the electrical conductivity raised enormously compared to 
the low doping case as expected. As a consequence, the powerfactor rises absolutely 
but unfortunately no enhancement via strain was obtained. The strain-dependent 
behavior of the thermopower as shown in figure 4(e) compensates the behavior of the 
electrical conductivity. In the high doping-high temperature regime the thermopower 
shows a noticeable anisotropy between in- plane and cross-plane component. While the 
in-plane component S^^ confirms our analytical predictions for high doping (see black 
solid line in figure 5) and even shows the shifted maximum to compressive strain. 
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Figure 6. (color online) Anisotropic thermoelectric transport properties for fixed 
temperature and hole doping concentrations in dependence on compressive and 
tensile strain. Left panels (a)-(c)) correspond to an hole doping of 1.6 ■ lO"'^ h/uc 
at a temperature of 100 K, while right panels refer to an hole doping of 0.16 h/uc at 
a temperature of 900 K. On the left axis of each figure the relative value compared 
to the unstrained case is shown, while on the right axis absolute values are given. 
Electrical conductivity ((a),(d)) and powerfactor ((c), (f)) are presented in units 
of the relaxation time t. 



the cross-plane component S'^does not follow the analytical model. This might be 
explained by multiband effects and nonparabolic bands, with the latter being more 
relevant in the cross-plane direction. The overall resulting powerfactor summarized in 
figure 4(f) shows however no significant anisotropy. Furthermore no enhancement by 
biaxial strain could be obtained, in contrast a decrease of about 20% occurs. 



3.2.2. Hole doping We will now focus in more detail on the thermolectric properties 
of hole doped biaxially strained silicon as presented in figure 6. Under full relativistic 
treatment the three valence bands next to the VBM at the F-point are the heavy hole 
(HH), light hole (LH), and spin-orbit split-off (SO) hole. While the HH and LH are 
degenerate, the SO lies 44 meV apart (see figure 3 middle panel). Under biaxial tensile 
or compressive strain, the valence bands become highly anisotropic and a crossover 
between bands occurs so that they even loose their original meaning[46]. It was shown 
that mechanical deformation-induced changes in the band structure offers potential 
for significant enhancement of the hole mobility [47]. Nevertheless, a straightforward 
explanation of the ab initio calculations as done for the electron doped case is not any 
more possible. Actually tensile and compressive biaxial strain does not only cancel 
the degeneracy of the heavy and light hole bands, which will cause reduced intervalley 
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phonon scattering, furthermore it also leads to a smaller effective conductivity mass 
and a further depletion of the uppermost hole band [14, 48]. Similarly to the electron 
doped case in figure 6 the thermoelectric transport properties for hole doped silicon 
under the influence of biaxial strain are shown for a fixed low doping, low temperature 
in figure 6(a)-(c) and fixed high doping, high temperature regime in figure 6(d)-(f). 
At low temperature and slight doping an increase of the electrical conductivity was 
found for tensile as well as for compressive strain for the in-plane component a^^ and 
the cross-plane component a ^ , while biaxial compressive strain tends to favor the 
enhancement of a. As shown in figure 6(b) the thermopower for hole-doped silicon 
again experiences a drop of nearly 7% under compressive and tensile strain. A possible 
explanation for the almost symmetric drop of the thermopower under compressive and 
tensile strain might be again related to the number of bands being occupied. At small 
doping for compressive and for tensile biaxial strain the primarily occupied HH and 
LH split and only the upper hole band is depleted and dominates the character of 
the transport properties [14, 48]. Extending our analytical findings for the electron 
doped case one would expect, that -Stcns. /compress. /so = (y)"^^"^ = 0.63. Through the 
counteracting behavior of electrical conductivity and thermopower under strain, again 
no enhancement of the powerfactor could be found (cf. figure 6(c)). Only under strong 
compressive biaxial strain a significant enhancement is visible in the low temperature/ 
low doping case. In the high doping and temperature regime not only the upper 
hole band plays an important role in transport, furthermore the former HH, LH and 
SO have to be mentioned. As shown in figure 6(d) an enhancement of around 10% 
compared to the unstrained case can be found for the cross-plane component under 
small tensile strain. For the in-plane electrical conductivity a^^ under compressive strain 
only a marginal influence on the strain can be reported. The thermopower shows a 
small anisotropy of the in-plane and cross-plane components. While the thermopower 
is mainly decreased for compressive or tensile strain, we again see a broadening of the 
thermopower drop in dependence on the applied strain with respect to the low doping 
regime. Beside the absolute values of the powerfactor being around 30% smaller than 
in the electron doped case (comp. figure 4(f) and 6(f)), a slight enhancement under 
thermoelectrically relevant doping and temperature conditions could be found for the 
cross-plane powerfactor PF^under small tensile strain (see figure 6(f)). 

3. 3. Optimization of the powerfactor by doping 

While up to this point the powerfactor and the incorporated thermoelectric transport 
properties were studied for fixed temperature and charge carrier concentration in 
dependence on the applied biaxial strain, we want to gain further insight in the 
doping dependence. Therefore the amount of charge carrier concentration to optimize 
the powerfactor PF||and PF^at given strain and fixed temperature of 900K was 
determined. This temperature seems to be a common temperature for thermoelectric 
application of silicon based devices. Figures 7 and 8 represent the results for the 
electron and hole doped case, respectively. 

From figure 7 it becomes clear, that an enhancement of the powerfactor by 
sufficiently high electron doping can't be stated. With increasing biaxial strain the 
in-plane and cross-plane thermopower decreases compared to the unstrained case. It 
is worth mentioning, that the charge carrier concentration has to be increased for 
increasing tensile and compressive biaxial strain to achieve the optimal powerfactor 
under the certain strain condition. Nevertheless, even for a raised optimal charge 
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Figure 7. (color online) Anisotropic power factor optimized by the carrier 
concentration at given strain state for electron-doped silicon. The in-plane (cross- 
plane) powerfactors are drawn as red circles (blue diamonds). Lines are shown to 
guide the eyes. The temperature is fixed at 900 K. 




Figure 8. (color online) Same as figure 7 but optimized for hole doping. 



carrier concentration the powerfactor drops to about 80% of the value of the unstrained 
case for the largest strain values considered here. As an interesting fact one can 
see that under tensile strain the cross-plane powerfactor PF^is always larger than 
the in-plane component PF|| , while under compressive strain it is the other way 
around and PFj^is smaller than PF||. In contrast to heavily electron-doped silicon 
an enhancement of the powerfactor could be found for hole doped silicon, as shown in 
figure 8. Nevertheless, the enhancement is limited to the cross-plane contribution PF^ . 
Here we find an enhanced PFj^undcr small tensile strain of ^°-/ao = 0.2%... 0.4%, 
while for tensile strain of 0.6% a value of the cross-plane powerfactor similar to the 
unstrained case is reached. We note, that the charge carrier concentrations, which are 
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necessary to optimize the powerfactor in the hole doping case are about three times 
smaher than in the related electron doping case. Even though an enhancement of the 
cross-plane powerfactor under optimized hole doping and tensile strain can be found, 
unfortunately the absolute values of the powerfactor are sufficiently smaller than the 
absolute value of the PF||and PF^of electron-doped silicon under all strain condition 
examined here. 



4. Tendencies on figure of merit 



To evaluate tendencies of the figure of merit we are going to include and discuss 
experimental results for the lattice thermal conductivity Uph, which adds up to our 
calculated electronic thermal conductivity k^i for the total thermal conductivity k. 
Involving the transport distribution function in eq. 4 the electronic contribution to 
the thermal conductivity is calculated as 

^el ~ 7^{C' ' - ■ (12) 

For bulk silicon, porous silicon and thin films various temperature- and doping- 
dependent measurements are available [49, 50, 51, 12, 52]. It is well known, that 
the lattice thermal conductivity of silicon is strongly depends on temperature. While 
at room temperature a value of 87 ^/mK was reported, this value decreases to 
36 ^/mK at 900 K and 23W'/mK at 1200 K [12, 50] making silicon an high temperature 
thermoelectric. Even at high temperatures and applicable doping the electronic 
contribution to the thermal conductivity is only a few percent of the total thermal 
conductivity. For bulk silicon it is reported, that the heat conduction is impeded by 
higher dopant concentrations [51]. Here the lattice thermal conductivity reduction 
due to the scattering by dopant ions overbalances the increase of the electronic 
thermal conductivity. Due to this behaviour the total thermal conductivity k is weakly 
depends on doping for temperatures clearly above BOOK [49, 52]. Nevertheless, for 
nanostructured silicon, e.g. in silicon-based superlattices [11, 53] or nanoparticle bulk 
silicon [12], with a low lattice thermal conductivity of Kph < 3^/mK, it may occur 
that the electronic thermal conductivity contributes remarkably to the total thermal 
conductivity. At high doping levels above 1 • 10^" cm~^ and temperature of 900 K the 
electronic contribution k^i can be around 1 — 3 ^/mK and therefore approx. 25 — 50% of 
the total thermal conductivity. While it was recently shown [54, 11], that tensile strain 
could lead to a reduction of the thermal conductivity by up to 15% in bulk silicon as 
well as in silicon thin films we did not include this in our estimation. In figure 9(a) the 
doping dependent figure of merit for three silicon based systems is shown. Wc assumed 
doping independent experimental lattice thermal conductivities stated for high doping 
rates, to include the reduction of Kph by ion scattering. Values of Kph = 36 ^/mK [12], 
Kph = 3.8^/mK[12] and Kph = 1.9^/mK[ll] were used for single crystalline bulk 
silicon, nanostructured bulk silicon and a Si/Ge superlattice, respectively. Here we 
consider only electron doping, which appeared to be most promising. The calculated 
doping dependent electronic part k^i is shown in the inset of fig. 9 for three different 
temperatures. For bulk silicon (solid line in fig. 9) we find a broad maximum for the 
figure of merit ZT « 0.1 at 900K and high doping levels of iV « 6.5 . . . 18 • 10^" cm"'', 
which is in good agreement to experiments by Bux et al. [12]. Compared to the related 
powerfactor (see fig. 1) the maximum shifts to lower doping concentrations due to the 
linear increase of Kei at higher doping rates, and therefore a stronger decrease of 
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Figure 9. Calculated doping-dependent figure of merit for n-doped bulk silicon 
at 900K (solid line), nano-bulk silicon at 900K (dashed line) and strained Si/Ge 
superlattice at 300K (dotted line). The inset shows the doping-dependent 
electronic contribution to the total thermal conductivity for 300K, 500K and 
900K. Experimental data (open square and circle) from rcf. [12] arc given for 
comparison. Note the logarithmic axes. 



ZT at high doping rates. If a reduction of the lattice thermal conductivity can be 
achieved by nano inclusions as reported in ref. [12] a remarkable ZT of about 0.5 
at A'' = 3 • 10^° cm~^ at 900K could be obtained making silicon an interesting high 
temperature thermoelectric (dashed line in fig. 9). Nevertheless, at room temperature 
the lowest thermal conductivities have been stated for silicon and germanium based 
superlattices [11, 53]. Therefore the dotted line fig. 9 estimates the figure of merit for 
a Si/Ge superlattice with a period of 15 nm assuming no degradation of the electronic 
transport by the hetcrostructure. Here, at 300K, the maxium figure of merit is found 
around doping levels of 7 • lO^^cm"'^ with values of about 0.1, which are comparable 
to single crystalline silicon at 900K at ten times larger doping levels. 

5. Conclusion 

In conclusion, the thermoelectric transport properties of biaxially strained silicon were 
studied in detail with respect to a possible enhancement of the powerfactor. Two 
different doping and temperature regimes were analyzed in detail. A low doping 
(TV « 1 ■ 10^^ cm""^) and low temperature regime (T=100K) suitable for metal-oxide- 
semiconductor device applications and a heavy doping (TV « 1 ■ 10^°cm~^) and high 
temperature regime (T=900K) suitable for silicon based thermoelectric modules. 

It was shown tliat the electronic transport properties, namely the electrical 
conductivity a and the thermopower S are highly sensitive against strain. Nevertheless 
it was found that strain-induced effects in a and S compensate each other and no 
remarkable enhancement of the powerfactor by either compressive or tensile biaxial 
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strain can be reached. On the other hand a reduction of the powerfactor of upto 20% 
(30%) under electron (hole) doping due to the influence of strain was found. The latter 
was assigned to a large extent to band-structure related redistribution of electrons. 
Estimations for the figure of merit under electron doping arc given. 

As a general result we showed that when degenerate VBM or CBM exist the 
thermopower decreases with deformation due to a redistribution of electrons in 
energetically lifted valleys. This could be the general explanation of the reduction 
of the electronic thermoelectric properties in strain-influenced hcterostructures. 
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